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Abstract 

Consider N = n\ + ri2 + ■ ■ ■ + n p non-intersecting Brownian motions on 
the real line, starting from the origin at t = 0, with rii particles forced to 
reach p distinct target points at time t = 1, with f3\ < /?2 < • • • < (3 P - 
This can be viewed as a diffusion process in a sector of R . This work shows 
q that the transition probability, that is the probability for the particles to 

pass through windows Ej~ at times tk, satisfies, in a new set of variables, 
a non- linear PDE which can be expressed as a near-Wronskian; that is a 
determinant of a matrix of size p + 1 , with each row being a derivative of 
the previous, except for the last column. It is an interesting open question 
to understand those equations from a more probabilistic point of view. 

As an application of these equations, let the number of particles forced 
to the extreme points (3\ and (3 P tend to infinity; keep the number of par- 
ticles forced to intermediate points fixed (inliers), but let the target points 
themselves go to infinity according to a proper scale. A new critical process 
appears at the point of bifurcation, where the bulk of the particles forced 
to —\/n depart from those going to ^fn. These statistical fluctuations near 
^vq that point of bifurcation are specified by a kernel, which is a rational per- 

K~) turbation of the Pearcey kernel. This work also shows that such equations 

are an essential tool in obtaining certain asymptotic results. Finally, the 
paper contains a conjecture. 
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1 Introduction 

Consider N non-intersecting Brownian motions X\(t) < ^(i) < • • • < x N{t) on R 
(Dyson's Brownian motions), all starting at source points 71 < 7 2 < ■ ■ ■ < 7^ at 
time t = and forced to target points S± < 82 < ■ ■ ■ < at t — 1. According 



to the Karlin-McGregor formula p2], the probability that the N particles pass 
through the subsets Ei, E 2 , . . . , E m C R respectively at times < ti < t 2 < 
• • • < t m < 1 is given by (setting t := and t m +i '■= 1)? 

/ m 

P p| {all Xi(t k ) G £ fc } 
\fc=i 

1 „ jv . at r N 

= 5^ / II ^ / I I • • • / I I det (P(^ - *05 7i, «f } ))k^ 

x det(p(t 2 - *il «i , 4 ))i<ij<iv ■ ■ ■ det(p(t m+1 - t m ; u[ m \ 5 i )) 1<iji<iV (1-1) 



for j 



N 



2 



where p(t, x, y) denotes the standard Brownian transition probability, 



P(t,x,y) := —= e * . (1.2) 



There has been a great deal of interest in non-intersecting Brownian motions 
and especially in some critical infinite-dimensional diffusions arising when the 
number of particles iV — > oo. This in turn has been motivated by random matrix 
theory and Dyson's observation [13] that letting the entries of GUE matrices 
run according to independent Ornstein-Uhlenbeck processes leads to such non- 
intersecting Brownian motions for the random eigenvalues of the matrix. 



When some source points and some target points coincide, the formula (1.1) 
for the probability must be adapted by taking appropriate limits; see [T7J HEJ [TUJ, 
[7J. In this paper, we consider the situation where the source points all coincide 
with 0, while some target points may coincide. Consider thus N = ni+n 2 +- ■ ■+n p 
non-intersecting Brownian motions starting from the origin at t — 0, with 
particles forced to reach p distinct target points (3i at time t = 1, with fix < (3 2 < 
• ■ • < (3 P in R; see Figure [TJ 

Given positive integers n = (m, . . . , n p ), given m subsets Ex, ... , E m C M and 
times t = < t\ < t 2 < ■ ■ ■ < tm < t m+1 = 1, this paper deals with the 
probability as in ( |l.3[ ) below (i.e., the probability for the particles to 

pass through the windows Ek at times tk)', as is well-known, (see [201 HQl EES E]), 
F^\t,K) can also be viewed as the probability for the eigenvalues of a chain of 
m coupled Hermitian random matrices, after some change of variables: 

/ m all ar f (0) = 0; 

t, E) := P I < all Xi(tk) E Ek \ n j paths end up at (3j at t = 1 



n \ i 



k=i for 1 ^ j ^ p 

Y[dM k 



. m m — 1 

If -§tr E M*-2 E c k M k M k+1 -2AM, 

\k=l k=l 



Z n Jspec(M fe )eE fc 



2 

e 



k=i 



=: P^(c,E). (1.3) 
The change of variables is given by the following formulae 



A := diag(&i, . . . , bi,'b 2 , . . . , b 2 , . . . , 'b p , . . . , b p ), with b e = x \ , 2 ^"\, / m - r fa 

\L — i rn )\i-—i m -i) 



5 E = Ex x . . . x E m . 

6 For m = 1, the matrix integral above becomes a one-matrix integral with external potential. 

it 



The change of variables below becomes: be = J j~tPti E = E.i 



'(i-t)' 
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^ ^ ; 2(t fc+1 — t k _i) ^ (t fc+2 — t k+ i)(t k — tfc_i) 

for £ = 1, . . . , p and k — 1, . . . ,m. It is quite natural to impose a linear constraint 
on the rescaled target points j3i, . . . ,/3 p , namely 

v v 

j^Ktfii = 0, with 2_^ K i = 1, set k := — 1. (1.5) 

Of course, the same relation holds for the Vs. For instance, a typical situation is 
to take (3\ = —f3 p and have all the remaining target points in arbitrary position 
between f3\ and (3 P . This case will be discussed in Section [8j 

The natural initial or rather "final condition 11 for the transition probability 
(1.3) is given by what happens when t m — > 1, keeping ti, . . . , i m -i, away from 
or 1; namely, 

lim Pf E) = 0, when £ m £{/V ••,&}■ (1-6) 

It is also known (see ([E])) that the probability above P« (ti, ■ ■ ■ , t m , EiX. . . X 
E m ) = det(l — X_ c (x)H^\x,y)X^ c (y)) can be expressed as a matrix Fredholm 

determinant of a matrix kernel H] 



27rV(l-*fc)(l-*£) 

'U-/3,, 



t_kXl,2xV_ 
l-t h ~<~l-tu 



><n 



V-p r J U-V 



7 The Fredholm determinant of a matrix kernel f/^t . (x, y) := XE t {x)H ti t j (x, y)xE j (y)' 
det (l-z(H ut M^r n ) 



„=1 0<r 4 <n ■"<■ 1 1 \ 1<3<^ / 1<fe ^ r 

where the n-fold integral in each term above is taken over the range 



-00 < a! 1 ' ^ . . . ^ o$ < 00 



n 



(m) ~ ~ (rn) 

1 m j 



4 




t-4. 



Figure 1: Non-intersecting Brownian motions 



0, for t k > 



, dy e e 1 "** for £ fc < ^ 

(1.7) 

where C is a closed contour enclosing all the points /3 r , which is to the left of the 
line := L + z'R by picking L large enough, guaranteeing 3F?e(C7 — V) > 0. 

These non- intersecting Brownian motions Xi(t) < ... < £tv(£) describe a 
diffusion process in a sector {xi < X2 < ■ ■ ■ < xn} of M. N and thus satisfy a 
diffusion equation. When the number N of particles tends to oo, the transition 
probability would have to satisfy an " infinite- dimensional diffusion equation", 
which however would be very difficult to use. The main result of this paper is to 
show that this transition probability P„(c, E) satisfies a non-linear PDE in the 
boundary points of Ei, . . . , E m , the target points bi, . . . , b p , and the couplings 
Ci, . . . , c m -i- It is the determinant of a certain matrix of size p + 1; p being the 
number of target points; so, when the number of particles tends to oo, the form 
of this equation remains the same, which will be exploited in the limit discussed 



in Theorem 1.3 Moreover, this determinant misses to be a Wronskian by the last 
column only. 

The PDE for the transition probability stems largely from integrable theory; 
this at least is our approach in the present paper. The integrable theory behind 
non-intersecting Brownian motions has been developed by us in j8] ; the latter con- 
tains many different ingredients; among them, multi-component KP hierarchies 
E] and multiple-orthogonal polynomials jU [HI HO] ■ It is - in our opinion - an 
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interesting open question to understand the PDE from a more probabilistic point 
of view and to use more conventional probabilistic tools to derive them. 

Throughout the paper, we shall use, without further warning, the following 
notation: (i) The inverse of the following Jacobi matrix will play an important 
role: 



/-I 



3 ■ -- 



Cl 



\ 







Cl 



-1 



Cm— 1 







V 



(1.8) 



c m -i -1 / 

(ii) For any given vector u — (u±, . . . , u a ), we denote by 



i=l 



d_ 



i=i 



cW 



(1.9) 



In particular, given any interval or disjoint union of intervals E = UjLj^i-i, ^2i], 
we denote by 



d E :~- 
e E :-- 



V 2r 8 

2si=l dzi 



sum of partials in the 
boundary points of E 
Euler operator in the \ _ <> 
boundary points of E ' ^«=i * dzi • 



;i.io) 



(iii) In view of the Theorem below, given b = (bi, . . . , b p -i) and subsets E iy define 
the linear differential operators: 



dt 

£b 
BO 



p-1 



flT:=0. implying £4" = 0, 



db 



e=i 



i=i 



X 



Jii for£ = 0, 
J" mi for 1 < £ ^ p, 



d 



£ El o~i,m£ b Ci , e m :— £ Em e 6 c m _ r 
For brevity in the statement of the Theorem, set ' :— d — Y^v 3ud E . . 
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(i.ii) 



Theorem 1.1 The probability F n : = P^(c,E), as in with the linear con- 

straint (1.5) on the rescaled target points, satisfies a non-linear PDE in the bound- 
ary points of the subsets Ex, ... , E m and in the target points b\, ... , b p ; it is given 
by the determinant of a (p + 1) x (p + 1) matrix, nearly a Wronskian for the 
operator ' := d , 



( 



det 





F 2 


F 3 ■ 


■ F p 


F{ 


Fi 




■ f; 


Fi' 


F» 


FS ■ 


. F" 

p 


f (p) 




F (p) 


F (P) 



G \ 

G 2 

G P J 



;i.i2) 



where the Fi and Ge are given by 
■■ -d d e \n¥ n - 



F t 

Ge+i 



i 



d o G e +j2(d o yFAd 



H 



a) 



-a 



H 



(2)' 



G t , 



0. 



i=l 



H, 



(i) 

t 

(2) 



(Ki(6 1>m - e m )d + 2J lm df) lnP n + C e 
(Si, m - + 2J lm bid )d e lnP n , 



;i.i3) 



with 



bp — b; 



(1.14) 



The final condition (1.6) translates into an " initial condition" near c m _i — > 
and 6^ — > oo, wpon ■usmg t/ie /act £/ta£ 



Cm-l - VI - t m , Cm-ibt ~ 0(1). 



As a special case, we consider the one-time probability ¥^(t, E) for < i 



t < 1. For this case, (1.3) becomes a one-matrix model with external potential 
P n := F^(E), thus with no coupling. The exp ressio ns for (1.13) can be replaced 
by simpler expressions; note that the ifj in (1.15) below are not obtained from 



the as in (1 1 . 13h , by setting m — 1; in fact, a further simplification occurs in 
the equations; also the functions Ge are only specializations of the above Ge up 
to a sign — (— 1)^ and ' now denotes Be instead of <% = —8e- In this statement, 



in accord with (1.9): 



we use the operator df^ as in (1.11), and we use the following simple operator, 



Ee — £bi with Ef, 



p-i 
i 



d_ 
dbi 
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Corollary 1.2 When m = 1 (th e one-time case), then lnP n = \nF^(E) satisfies 
the same non-linear PDE (1.12), but with simpler expressions Fi and and 



with ' = d E , 

F t : = (di e) + K e d^jd E \nF n + n e , 

Hf ] := (-K l d E e+{K l {e-l) + 2){dl l) + ^dj) lnP n + Q, 

: (l-e + 2b e d E ) (di e) + K t d^]n¥ n , (1.15) 



H. 



(2) 



G e+1 := d^ + J^idjFiid^-dP^-), G = 



i=l 




Ct := —2ni 



In section [7j we shall work out two examples, immediate applications of the 
equations in Theorem |1.1| and Corollary 1.2| In the first example, we describe 



nonintersecting Brownian motions, leaving from and forced back to 0. The 
second example deals with the situation of several target points with the extreme 
ones being symmetric with regard to the origin. That model will also be used 
later in Section [HI 

Pearcey process with inliers: In section |8j we consider non- intersecting Brow- 
nian motions leaving from and forced to p target points at time t = 1, with 
the only condition that the left-most and right-most target points are symmetric 
with respect to the origin, with p — 2 intermediate target points thrown in totally 
arbitrarily; it is convenient to rename the target points (3\ < . . . < (3 P , as follows: 

a < — £i < . . . < — c p _2 < — a 

(1.16) 

n + n x ... n p _ 2 n_ 

with the corresponding number of particles forced to those points at time t — 1. 
The purpose of this section is to identify the critical process obtained by letting 
n := n + = n_ — > oo and by rescaling a and the c, accordingly, while keeping 
Tlx, ... , n P -2 fixed. We let a go to — oo like —^fn and —a to oo like y/n. The target 
points — Ci, . . . , — 5p-2 of the inliers move to oo as well, but at a much slower rate, 

namely like — ue (l) 1 4 . A new process will appear at the point of bifurcation, 
where the bulk of the particles forced to —y/n depart from those going to \fn, 
namely the Pearcey process with inliers, which generalizes the Pearcey process 
found by C. Tracy and H. Widom [TS]. It describes the statistical fluctuations 
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near that point of bifurcation; it will be sensitive to the presence of inliers and 
will be different in the absence of inliers (Pearcey process). We will compute the 
kernel governing the transition probabilities and also apply the formulae obtained 
in Corollary 1.2 to compute a PDE for the gap probability, which, to our surprise, 



appears to be an exact p x p Wronskian. This is the content of Theorem 1.3 



Theorem 1.3 Pick times T\ < ... < T m , subsets Ej C R for j = l,...,m 
and parameters U£ for i = 1, ... ,p — 2. Consider 2n + X]?=i n £ non-intersecting 
Brownian motions, such that 

(i) all particles leave from at time t = 0, 

(ii) n = n± particles are forced to ±\/n at time t = 1, 

(Hi) ne paths are forced to poinh 8 —ui (f) 1 ^ 4 at time t = 1 (1 ^ I ^ p). 

Then the following Brownian motion limit holds for the gap probability, about 
time t = 1/2, keeping ni fixed, 

1 TV \ E c - 




all Xi(- H 1 -i==] 

V2 4x/W 



4(n/2)V 



i 



F v ( U1 ,...,u p _ 2 ) ^Q{p( Ti ) n Ej = 0} 



<l>3 



dct ( 1 — yX E K^_ T X B _ ) ), (1.17) 



where this probability is given by the Fredholm determinant of the Pearcey matrix 
kernel with inliers, which is a rational perturbation of the customary Pearcey 
kerne¥\ namely 



K s,t( x , Y ; ui,...,up- 2 ) 
1 




;i.i8) 



8 Note that those points belong to the interval [— y/n, y/n\ for large enough n. 
9 X stands for the contour 

^ .0* 
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The log of the gap probability (E = E\ x • • • x E n 



,...,r m ; Ul ,..., u p _ 2 ; E) := lnP^ («i.-..,«p-») ^ n E . = 0} 

satisfies a partial differential equation, which is a p x p Wronskian with respect 
to the operator d E = Y1T=\ &e '■ 



d 2 — 
' ' ' ' ' E du 



p-2 



0. 



1.19) 



when 10 



X := - e u + 2e T - 2)<f Q + 4d u %d T Q + 8d 3 T Q - Hc\d T Q + 4 {c\d T Q, % 



For one-time (m = 1), the expression X reads as follows: 
X := (e E -e u - 2r^ - 2)£Q + Ad u d E + 8^ + 4 # 



1.20) 



• (1-21) 



Remark: The term e„d 2 Q could be omitted in the definition of X, since it is 



a linear combination of (p— 1) columns in the matrix (1.19) (from the second Xu± 
to the (p— l)st column x« p _ 2 ). We nevertheless keep this term in the expression, 



in view of Conjecture 1.5 



In the absence of inliers, one obtains, in particular, the PDE for the transition 



probability of the Pearcey process: it is a 2 x 2 Wronskian with X as in (1.20) 



and (1.21), but without the w-partials. In [3], it is shown that the transition 
probability of the Pearcey process satisfies the simpler equation X = 0. 

Corollary 1.4 In the absence of inliers (p = 2), 

Q{P(r,) n E d = 0} 



satisfies 



+ 2e T - 2)d 2 Q + 8«9 T 3 Q - Ad K d E d T Q + 4 {c^Q, d t 



10 



Remember for u = {u\, . . . , Up-2) and r = (ri,...,T m ), one has d u = 



% ■= Ll T iO E 



d T = j^r-One also needs Ee '■= Ei™ £ £i an d the mixed time-space derivative 
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-vflf) 








Figure 2: Pearcey process with inliers 



and for the one-time case (m = l), 





0. 



We now formulate a conjecture, stating that, even with inliers, the equation 



for the transition probability reads X = 0, where X is given by (1.20) and (1.21 ): 



Conjecture 1.5 Even with inliers (p > 2), we conjecture that the function 

tm N 
Q{P(r J )n J E J = 0} 

satisfies 

X = (e E -e u + 2e T ~2)d^Q+Ad u cld T Q + 8dfQ-4dAdrQ + ^{d E d T Q,d E ^ ) () 

and for the one-time case (m = l), 




2^ _n 



1.22) 



X= (^-^-2^-2)^+4^^+8^+41^^,^. 



0. (1.23) 



The PDE's play a prominent role in obtaining certain approximations which 
would be very hard to obtain without that technology. An example will be given 
here, without proof, for the Pearcey process without inliers. At the point of 
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bifurcation, mentioned above, there appears a cusp in the Pearcey scale £ = 
±|^(3t) 3 / 2 , such that, roughly speaking, most Pearcey process paths stay com- 
pletely to the left or to the right of this cusp. Upon comparing the Pearcey 
process with, say, the right branch of the cusp in the new (crude) space-scale 
(3-r) 1 / 6 , and letting two different times r% and r 2 tend to oo in a very specific way, 
one is led to the so-called Airy process A(t). The exact approximation is given 
in the Theorem below taken from pQ: 



Theorem 1.6 Leiri,r 2 

t~2 - n 

2(t 2 - £1) 



oo ; such that 
1/3 , h — ti 



(3n)^ + 



(3n)V3 3n 



this specifies two new times t% and t 2 . The following approximation, far out along 
the cusp, of the Pearcey process by the Airy process holds: 




2_ 

27 



(3r, 



A 3/2 



P( [){A(ti)ru 



n (-Ei 



-E, 



\i=l 




Remark: The O (r x 4 ^ 3 )-approximation, obtained via the PDE is much better than 
any rough estimate one might predict. Also one expects that, in this precise limit, 
the Pearcey process with inliers tends to the Airy process with outliers; see [2]. 



2 Non-intersecting Brownian motions and a chain 
of Coupled Random Matrices 

Setting 

1 1 1 



i~fc := t k+ i - t k and 



0~k tk — tk-l tk+l—tk 



for 1 < k < m, 



and taking in (1.1) the limit 7, — > 0, for i — 1, . . . , N, leads to 



/ 111 

P( p| {all Xi(t k )eE k } 

rn 

Ajv(Mi) ] [ 



\k=l 



x j (0) = 0, x j {l) = 6 j , 
forj = l,...,iV 



ZL .ie n 



k=l 



2u k;i u k+l;j 

det I e T k 



U k ; 1 



] J e "k du^i 

l<i,i<JV \<{<n 



(2.1) 
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where Ajv(wi) stands for the Vandermonde determinant in the variables ui = 
(ui-i, . . . , Ui-n). Notice that each of the sets of variables u±, . . . ,u m appears in 
exactly two of the determinants in the above integrand and that the other factors 
are insensitive to a permutation, for fixed k with 1 ^ k ^ m, of the variables 
Uk = ■ ■ • , Uk-N- Therefore, taking the limit u m+ i ;i = Si — > fa, for i = 1, . . . , N, 
with n e of the Si going to fa, namely w m +i ; i, . . . , u m+ i- ni — > fa, and so on, making 
m synchronized changes of variables, and using the symmetry of the integration 
ranges vis-a-vis these variables Uk-i, ■ ■ ■ , Uk-N, 



P 



m 

p| jail Xi(t k ) G E h 
k=i 



x j {0) = 0, (j = l,...,N), \ 
Xi(l) = • • • = x ni (l) = fa, 



XN-n p +l(l) 



x N (l) = (3 P J 



1 



1 
7/// 

1 



3 <V 



Aiv^i) n a 



oo 2 



i ~ (0 CO 



-+ £ 



1=1 



nil 



M«i)nK(*n 



i (<0 2 , "V 1 («) 00 , l 00 
> ^v, ' + > etui -v, i , . 



e fc=1 



„ m 

I In(v) TT 
J ^ N fc=l 



|tr( £ Af 2 -2 £ c fc M fe M fc+1 -2AM 71 



spec(M k )eE k 



fe=l 



I { du 



k;ii 



][ dv k;i , 



(2.2) 
(2.3) 



where the diagonal matrix A, c*., 6^ and were defined in (1.4) or alternatively 
expressed below in terms of the a^s and Tfc's. The last integration is taken over 
Hermitian matrices, with spec(M&) e Also the change of integration variables 
u ah l— ^jh above is given by 




For 
by 



m and for £ 



-, oe — pi, hik 



Tk 



T,, 




, p, the vector u 



(i) 

k 



,(0 



<?k 



l u M>---> u k;n £ J is defined 



(Wfe; 



1) 



(1) (1) (2) (2) (p) (p) >, 

• • • ' K fc;ni' U fc;H • • • > M fc;n2 ' ' • ■ ' U k;\i ' • • ' U k;n v ) 
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Concerning the Jacobi matrix (1.8), one needs the following formulas for deriva- 
tives of J] they can be shown by recurrence: 







^1 n >smm 
OCi 



yJrr 



-2 r 2 







-TO— 1 



dc 



J, 



ml 



(2.4) 



m— 1 



3 Integrable deformations 



In this section, we introduce a time deformation I n (v) of the integrand I n (v), 



introduced in (2.3). The deformation is chosen such that the resulting integral 



is on the one hand a solution to the multi-component KP hierarchy (see [8] 



and Proposition 3.1 below) and satisfies on the other hand a set of Virasoro 
constraints. We will impose on the rescaled target points b±, .. . ,b p , which we 
henceforth denote by . . . , vf \ a non-trivial linear constraint 



0. 



(3.1) 



Without loss of generality, we may assume (upon reordering) that k p ^ and 



impose if Ei K t that Yl\ 



1; also define k := —1. Thus, the non- 



deformed integral which we will consider is 



3 Ar 



wIew^oII dVk 

k=l 



(3.2) 



The integrand I n (v) will be deformed by four sets of parameters: (i) A first set, 
denoted by , . . . , b ^ , deforms the parameters b^p 
same constraint (3.1) as the parameters b\ , namely 11 



2 



0. 



(3.3) 



(ii) A second set of deformations consists of parameters corresponding to the KP 
time variables; they are denoted by Sr ( r G ^>o) for the parameters going with 
the starting point of the Brownian motion and si (1 ^ £ ^ P and r G Z >0 ) 
for the parameters going with the £-th end point of the Brownian motion, (iii) 
There is furthermore a set of parameters 7^ (2 ^ k ^ m — 1 and r G Z>o) 



by Ei K i b il 



n The combination of the two constraints (3.1 1 and (3.3) will in the formulas below be denoted 
0. 
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going with the intermediate times t 2} . . . , £ m -i and (iv) a set of parameters c, 



(fc) 



r,q 



(k = 1, . . . , m — 1 anc£j (r, g) > (1, 1)), going with consecutive times t^, t k+1 . 

For n = (jix, . . . , n p ) and E = E\ x i? 2 x • • • x £7 m , where each ^ is the union 
of a finite number of intervals in R, define 



r n (E) := 



4W 



n ^ 



(3.4) 



fe=i 



where 



In{v) = I n {v) X 



with 

Uv) 



P ni h W„W 



«> E rfMf -«>/)+ E E 



nn 

£=1 8=1 



Cr 9 U fc;j W fc+l;j 



fc=l (r, 9 )>(l,l) 



M e &=i fc=i 



i=l 



We denote by £ the locus corresponding to setting all deformation parameters 
equal to zero, so that I r 



'r ; • • • ; °r 



0, r G Z 



>0) 



6 W b (p) _ o 

/r , • • • , 7r 



0, r G Z 



> . 



>0; 



.(1) 



n (m-l) 



-0, (r,g)>(l,l) J 
We list a number of operator identities, valid when acting on r n (E) 



(3.5) 



5 



96 



h 



£=1 ^ 
9 



P-1 



d Kf d 

+ -— y i^e^ P -i,h = i,2, 



(0 



56 



Me {1,2}, 



-(1-5. 



9 



96 



(0 



A d ^ d 



P Q 

1=1 



(3.6) 
(3.7) 



/i = 1,2, 1 < I *C p, (3.8) 



12 



The inequality (r, q) > (1,1) means by definition that r J? 1, <? ^ 1 and (r, q) ^ (1,1)- 
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where for h — 1,2 and 1 ^ £ ^ p we define 

p-i 



P-i 



°b h ■- 



:i-s t 



PJ db {i) +Ke ^db (i) 

uu h i=l uu h i=l 



d 



db 



implying 



(3.9) 



(3.10) 



t=l 



Using Yfi= i K $h^ — OjOne first establishes identity (|3.6[) and then J3/7J), while 



the first equality in (3.8 



using Yj£=i — ^ an< i t 



is obtained by computing 5Ji=i ~~To fr° m (3-6) and by 

I 

ie identity (3.6). 



From section 7.3 in [S], it follows that r n (E) can be written as 
r n (E) = det 



/ ((yVl(y) I X j l/j(x))) 0|i<nj \ 



V ((y i( fp(y) I ^W)) ° 



5^ i < rip 
< j < JV 



(3.11) 



/ 



where 



exp ^-^ 2 + ^4°vj, 

exp f-^ 2 + ^+4V-E^V 

V 1 T>\ 



for 



,p, and where the inner product (• | •) is defined by 



(f(y) I g(x)) 



f(y)g(x)fi(x,y) dx dy, 



E 1 xE„ 



with 



fi(x,y) :-- 



m—l 



exp 

•■-> '-' k=2 



[-l w t+J2^ k)w i i >< 



r>l 



m—l 

exp 

k=l 



m—l 



c k w k w k+1 + 

(r,5)>(l,l) 



C rq W k W k+ 



i fj dw* 



fc=2 



tOi := a; and u; m := y. For m = 2 the latter formula for /x should be interpreted 
as fJ,(x,y) := 1, while fJ,(x,y) := <5(x — y)e x I 2 (the delta distribution) in the case 
of m — 1. 
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The above representation (3.11) of r n implies, in view of [El Prop. 6.2], that 
T n is a tau function of the p + 1 component KP hierarchy, in particular we have 
the following Proposition. 



Proposition 3.1 The function r n = r n (E), as in (3.4), satisfies for 1 ^ £ ^ p 

d 2 







ds ( ; 



(0) 



In 



T, 



n+e e 



lnr„ 







T, 



n—et 



d 2 



lnr„ 



ds\ 



it) 



In 



' n+eg 



ds^ds? 



lnr„. 



t. 



n—et 



d 2 



lnr„ 



(3.12) 



where n±e e = {n x , . . . ,n p ) ± e t := {n x , . . . ,ne-i,ne ± l,n e+1 , ...,n p ). 
Both equations will play an important role in Section [6] below. 



4 The Virasoro constraints 



Rem embering the definition (1.10) of the operators 8e and Ee and the definition 
p.9) of the operators df \ define for I — 1, . . . , p the operators: 



(o) 



p-i 



iTifc^E* — 2j7i m b 



(0 



9 



fc=i 



p-i 



i=l ^ 



2 



,fe=l 



(o) 



5 



(0 9 , 2 ^ fe(i) 9 



l(<) 



B. 



p-i 



p-i 



<9 



Cc m _i , =1 0&w i=1 ^ 



(4.1) 
(4.2) 
(4.3) 
.(4.4) 



We show in the following proposition, how the action of these operators on the 
tau function can be represented by time derivatives. 

Proposition 4.1 The integral r n (E), as in (3.4), satisfie^\ for I — 0, . . . ,p and 

h = l,2, 



Bf hxr n 







ds 



CO 



+ KtEjP I lnr n + k{T } 



h > 



(4.5) 



13 Recall that k = — 1 
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where 



(«) 



a = 0, 



1=1 

-Ji m Nsf ] - J mm Mb? ~ sf ] ) a ^ 0, 



(4.6) 



f E 



(a) 



riiUj 



iV(iV + l)/2 



??? 



a = and m > 1, 



(4.7) 



^^ni(rie + l)/2 a ^ and m > 1, 



and eac/i E^ a homogeneous first order differential operator in all deformation 
parameters, except for the deformation parameters b\ f , so that E^, = ; and 
moreover, for k = 1, 2 and /or £ = 1, . . . , p, 



(f) ' ^ 



5/1,2^1, 



<9 



(0) ' ^2 



<9 



(or 



(4- 



Proof: We give a detailed proof for the case of m = 2 (see remark 4.2 for the 



case of m > 2 and see remark 4.3 for the special case of m = 1). Then c m _i = C\. 

2 matrix 

1 A c N 



which we simply write as c. Also, J7 is the 2x2 matrix 

-l 

/ — i i- \ 

J 



-1 c 
c -1 



1 -c 2 Vc 1 



In this case, referring to (3.4), there are two sets of variables v\ and t>2, which we 

(k) 

denote by x and y, there are no deformation parameters 7r and there is a single 
set of deformation parameters c^ q , which we will denote by c rq . For E\,Ei C M., 
and taking into account the usual constraint Yle=i K lb\\ — 0, 



r n {E 1 ,E 2 ) :-- 



Yli=i n £ ] - JJe^xE. 



^■1 



I n (x,y) dxdy, 



where 



Tl£ 



I n (x,y) = A N (x)HlA m (y)H 

i=i \ t=i 



i a) 2 i w 2 . (i) a) , ,«) (i) 

e -2 x i ~5fi +ca: i »< + 6 i fi X 



(4.9) 



(4.10) 



\ w )+ e «v,*r J i/r 

(r,9)>(l,l) 



We first compute the action of the operators ds k and EE k on the tau function 
( 4.9[ ). We start with 8e 2 - Using the fundame ntal t heorem of calculus and the fact 
that J2f=i ^:^Jv(y) = 0, we compute from (4.10) that 



dE 2 T n 



N 



i=l y% 



E E -yf + <*? + tf> - E fc 4 £) (yf ) fe " 1 + 2 W 



(r,?)>(l,l) 



£=1^*1 f*! £=1 £=1 fc>2 Ub k-l 



E 



Crl" 



(9 



9 

k-1 
\ 



E gcr 



d 



1 Ub \ r^2 (r, 9 )>(i,i) r ' 9 1 



P 



/ : os i as i £=i 



(o) 1 / "' 

(»-, q 

q^t 2 

/' 

'""ft 



I n (x,y)dxdy 



a 0) +^ & 2 (£) 
1 fc>2 ""fc-l fcl UU 1 



Ed \ - (9 

r^2 OS r (r,<j)>(l,l) T ' q 

9^2 



where we have used the identity (3.7), which follows from the constraint Yli K iof\ = 
in the last step. The computation for 8e 1 is similar, but simpler: 

N 

E dxdy 



i=l 



EE 

i=\ 3 = 1 



+E*4 0) (*?)*- 1 



0.s 



9 ^ d 
- c E 



(0) 



=1 ds l 



, ^ 0) + E fc 4 0) — 



E rcrM?Y-\v? ) q )L(x,y)dxdy 



r, (0) Z^l^"^1r, (£) ' ^ 

fc>2 OS fc-l fcl ?>2 OS <? (r, «)>(!, 1) T ,q 



EE 



ci g - 



E 
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For the computation of the action of e E and e E on the tau function, note 



i=l 



t dxi / —' axi 

i=i i=i 



and so from (4.10), compute using (3.7) and the constraints Yli K ^i2 = 



2 



ff N d 
J J ^2o^_(y^n(x,y))dxdy 



f=i 



E E I -yf 2 + + *W - E + *Pvf 

e=i j=i \ fe^i 



Y, qc rq {xfy{yfy\ ]l n (x,y)dxdy 

(r,g)>(l,l) 



E 



Pa q P 

^ 9s W ' C ~dc 9s 

= 1 a6 2 f=l a *l 



n £ (n £ + l) , v ~ i) if J2 b W 9 ^ | J2J2 k "k ] " 



E 



V 



d 



4 (i) + E qCrqTj— \ In{x,y)dxdy 

1=1 ® s 2 (r,g)>(l,l) Cr? 



(4.11) 



p n/(n/ + i) <A a a . s 



E 



E 



E^dW+EE** 



p-1 



d ^ d 

db w 9Cr9 a^ 1 r " 

= 1 OT 2 (r,q)>(l,l) r? 



Similarly, 



N 



Ef xE? 



i=l 



9_ 



2J ■£ r (a; i / n (a;,2/)) ctecfy 



9 



rc 



aS fc (r,g)>(l,l) 
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In order to deduce (4.5) from th ese for mul as it suffices, for h — 1, to substitute 
in the first line the definitions (4.1), (4.2) for B-f and in the second line the 



expressions for d^ x T n and 9^ 2 r„ in 14 

'«o 



B\ 

A 



(0) N 



(0 



Ki 



J 



d El 



(I)) 



1=1 os l 

a 



k 
^ 



J 




\i=l 




^ + E<°>) 
«*(T« + E?) 



+ 



«oEi 



(0) N 



where we used (3.8) (for k — 1) in the third line, and where we set 



El' 



(0) d 
k 



ds 



(0) 



Ya=1 Sg^2 C l9lS) + S (»-,«) >(l,l) r C, 

Us n r > 9 



r^2 C rl g a (0) + £ 9) > (1. i) <?Crg 



r 9 9c r _x, g 
9 



9^2 



(9c 



r,g— 1 



(4-12) 

Thus we see that E^ and E^ are homogeneous first order differential operators 



^(0 



(i) 



M 

, . . . , a 1 , 



in the deformation parameters, and that they are independent of 
and of . . . , 63 , leading to the stated properties of E^ and sf' 1 . For k = 2, 
it suffices t o sub stitut e th e found expressions for ee 1 and Ee 2 , acting on r n , in the 
definitions Q| and O of ^ 0) and B? , to wit: 



o(0) 
i_> 2 1 n 











dc / 

N(N + 1) 



. ,(0) ^ 



«94 0) 



(r.5)>(l,l) 



<9 



rq dc. 



rq 



ds 



d N(N + V 
+ ^o^^ 



(0) 
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Recall that m — 2 and that kq = — 1. 
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and 



' a 



/ >A mini + 1) <A ^ (i) a ^ .9 



«=1 «=1 fc^l tyo fc (r,g)>(l,l) 

+ S H T n 



Ub 2 i=l 



wherd^l 



n (o) 1 ' ^'i Q c • 

fe^l aS *: (r,g)>(l,l) r? 

m d ^ d 



p 



i=l fc^l " S fc (r,g)>(l,l) rq 



Remark 4.2 For m > 2 the proof goes along the same line, but it has extra 
terms, coming from the deformation parameters 7^ . As it turns out, 

f) m ( ( p ~ l f) p W 

"7l i=l \ \ £=1 tf&i £=1 / / 

while before, mod 0(C), so the — ^yT n are only needed to solve 

for —70 r n in terms of the (Be, — 5, m (*))T n , but they do not enter into the actual 

ds\ 

solution of -^ij)T n mod 0(C). 



Remark 4.3 For m = 1 (one time) the proof of Proposition 4.1 is simpler, but a 



few adjustments are needed. Denoting the subset E 1 C K by setting k := — 1 



15 Notice that is independent of £ for 1 sC I ^ p. 
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and := d^' := 0, the operators B[ L> and B^' can for i — 0, . . . ,p, be written 



as 



i=l 
p-1 



(0 



(4.16) 



i=l 



while Ti. 



T< q) and S fc := £< a) are independent of a and take the simple form 



71 = JV S <°> + 5>(f>f - sf>> 



H= E 



5 



£=0 fc^2 fc— 1 



EE* s 

£=0 fc>l 



fe „ {(.)■ 



(4.17) 
(4.18) 



5 Virasoro constraints, restricted to the locus £ 



Restricting the operators Bi, T{ and Ej ((4.1) - (4.7)) to the locus C, defined by 



setting all deformation parameters equal to zero (see (3.5)), yields new operators 
for I — 0, . . . ,p, 



(0 



B\ 



o bl - n e 



i=i 
d 



x 



Ju for 







J mi for 1 < t < p 

p-i 



-£ Bl + Ci 



9ci 



<5l,m &l 



5 



m— 1 



<9c 



m—l 



(0 



for O 1. (5.1) 



while all E^ , defined in (|4_12j) and fl4.14p , restrict to zero, T 2 W = T 2 W for < £ ^ p 
and 



(0) 



J lm N(h), if = -JmmNih), for 1 *C £^ p, 



(5.2) 



where N(bi) := J^^i 71 ^! ■ ^ l ea ds, on the locus £, to the identities: 



Proposition 5.1 For £ = . . . ,p and h = 1,2, the following formulas hold on 
the locus C: 



In T n = B$ In r n - Kifj® 



ds 



h ' 



(5.3) 
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while for second derivatives and 
d 2 



,01., (Q lllT ™ 



ds\ ds\ 
d 2 

ds^ds? 

d 2 



B^B? \nr n + n e J lm 



lnr n = + 5 1 ^)B^\nr n -2J 2 m K,N{b l ) 



p, also on the locus C, 



(5.4) 



In r n = (B^ ] - Kt5i }m )&i' In r n - 2J lm (J mm K e N(b 1 ) + d^' In r n ). 



Proo f: T he first set of i dent ities (|5.3|) follows at once from restricting the identi- 
of Proposition 



ties (4.5 



4.1 



order differential operators. 



to the locus £ and using that -^py and Bjp are first 



The identities (5.4) involving second derivatives are 



shown as follows. Concerning the first one, observe from Proposition 4.1 that 



B$W ) In r n = B^BY-' In r n = B\ 
c c 



j(0) R W 



(0) 



d 



d 



(f) 



d 



+ K^ ) \Bf ) \nr n 



0sf 



da? 



lnr„ 







+ /t sf } ) lnr„ + k T^ 



(0) 



(0) 



d T-(°) 



In - 



J\m n ii 



where we used in the last equality the relations ^py^i — Jim n t (see (4.6)) and 



<L_ v (o) 



94' 

third one, we use that 



0, (see (4.8)). This yields the first identity in (5.4). To prove the 



i=l 



p-1 



6i ' 



as follows from (3.9), and 

dj h 



gW r (o) 



p-i 



5_ 

-KtJi m {2J mm + l)JV(6i) + KiJimN^x) = -2K t J lm J mm N^) x ), 



w<Cm _ 1 ^LJV(6 1 ) + K , Ji m V &« -^iV(6i; 
<9c m _i ^ ^ w 



by using (2.4), when m > 1, and B^T\ 
m — 1, so that 



-KeN(bi), by Remark 



4.3 



for 



B C0 T (0) 



-KeN(b 1 )(2J lm J mm - <5i 
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for all m. Using these identities, (4.1), (4.4), Proposition 4.1, (4.8) and (5.3), 
compute 



p-i 



BM 0) ^ r„ + 2J lm £ fl£> (&«) ° In r n 



i=l 



B. 



Co) 



2J lm d^ lnr, 



k S^ 0) lnr„ + K^(6i)(2J lmi 7 mm - <5i, m ) + 2J lm d^ lnr, 



(0) 



d 2 

o 2 



K(YSS > I lnr„ + /t^T, 



+ Ke N(b 1 )(2J lm J mm ~ S hm ) + 2J lm df} lnr ; 



lnr n 



In- 



+ K t N(bi)(2Ji m J mm - 5 1>m ) + 2Jx m d ( j® lnr„ 



In 



In - 



5(0) 



lnr n 



- JV(6i) I + 2,7i m (/s<JV(6i) J" mm + 9 6 ( f In r„ ) 



+ £i,m«^ u; lnr„ +2 l / lro (/t^iV(6i) l 7 mTO + ^7lnT n ), 

£ £ 1 C 



N{b x ){2J? m - S 1>r 



which yields the third relation (5.4). Using B 2 T\ 

I I c 

which follows from (4.3), (4.6) and(2.4), the second identity in ( |5.4 ) is proven in 

a similar fashion, using (4.8) and (5.3), namely 



B^B?\nr ri 



B. 



(o) 







K £ Y,f> I lnr,, 







ds? 




+ K t E?> ) B ( 2 0) lnr„| £ + K£ iV(6 1 )(2J 1 2 m - 5 1;t 



d 



dsf' \\ds 

d 2 



(0) 
2 



(o) 



+ «*JV(&i)(2J? m -<f 1 



dsi 'ds 



lnr„ 



9 



,(0) 



W ' 2 



lnr„ 



+ ^iV(6 1 )(2Ji 2 m -5 ljm ) 
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d 2 . 



ds^ ds 2 
d 2 



A 9 1 

Ol.m — OXTr, 



<9s 



lnr„ 



+ ^iV(6 1 )(2Ji 2 m -5 1 , m ) 



+ 2^iV(6 1 )Ji 2 m -,5 lim Sf ) lnr ) 



6 A PDE for the transition probability 



This section aims at proving Theorem |6.3[ which leads at once to Theorem 
In order to do so, we shall need two propositions: 

Proposition 6.1 For 1 ^ I ^ p, the function Xe := d^Bf lnr n satisfies the 
equation 



{Xe, Ft}g(o) 



9(0) 



TP 

tin .Ft 



(6.1) 



where (n = (ni, . . . , n p ) ) 



H, 



(i) 

i 

(2) 



c, 



lm j 



(K t (6i, m - e m )Bf ] + 2J lm d£>) lnP n + C t 



p-i 







£ E m ~ Cm-l 



do, 



m—l 







{£)■ 



0-n.T T fcW _ V" Ui 

\ i^e W — b 1 



(6.2) 
(6.3) 



Proof: From (5.3) and (5.2), one finds, along £, for 



■P- 



a4 0) 
9 



T, 



n—et 



In 



r, 



n+ei 



r, 



n—et 



(0) , T~n+e e 



B^> In 



r, 



n— e/ 



n-e t 



(6.4) 
(6.5) 
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A direct substitution of these formulas, as well as the formulas (5.3) and (5.4), in 



(3.12), leads, along £, for £ = 1, . . . ,p, to 

1 



K°) i^ Tn + e e 



B?> In 



i n-\-e£ 



'n-e e 



2j7lm&l 



+ 2fi^j7 mm ,^l 



(£?f + 8 ljm )B { ^ lnr n - 2J? m K t N{b 



(0) 



-2J lm (JmmKeN(h) + dif lnr, 



where F e := lnr n 



tions the term which contains 

m 



— UiJ\ m (see (5.4)). Eliminating from these equa- 



te) 

which can be done by applying B\ to the 



first equation and B\ to the second equation, and using that these operators 
commute, we get the single equation 



. . 2J lm bfF t - (B^> + 5 hm )BY> lnr n + 2J? m K t N(b lt 
W I c 



Ft 



{B? - K t 8 x>m )B^ lnr n - 2J lm (J mm K t N(b 1 ) + d^' lnr n 



Bl 



Using the fact that the derivative of a ratio amounts to a Wronskian, b y cle aring 



the denominator, and writing B^ as B^ = + KeS m (see (6.2) and (5.1)) and 



using the formula for Fe, one can rewrite the latter equation as 



(tff + 6 ltm + 2J lm b\ l >BDB\ l> lnr n + 2J* m (n t b? - «*JV(&i)), F t 



■Mn(.°hn(Q 



(6.6) 



^i(£m ~ $i,m)Bi'* ~ 2J\ m d^) lnr n — 2KiJi m J mm N(bi) 1 Fi 



g(0) 



Finally the integral r n (as in (3.4)), but integrated over the full range K, equals 
(see the Appendix) 



g n (c)e 2 U=i^ b i 



«) 2 



n (^-^ 



(6.7) 
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with g n (c) a function, depending on ci,...,c m _i and n only. Thus one has, 
restricted to £, 



lnr n (E) 



lnP n (E) + lnr„ 
p 



(6.8) 



J2Mb{ E) ) 2 + n inj \n(b[ j) -b?) + \ng(c). 



= 1 



When (6.8) is substituted in (6.6), a few terms will appear where lnr n (IR m ) is 



acted upon by a differential operator. We derive the formulas which will be used. 



First, it is clear that V n (R m ) 



0. Therefore, since 



b[°\b? 



0. 



F e = -Bf ) B[ e) lnr n (E 



lm 



BfB? \nF n (E) - n e J lm . 



Also, using d^bf* = m — 8^ valid for i — 1, . . . ,p, one computes 



^lnr„ 



b? - bf 



and therefore, since Bf ] lnr n (M m ) = 0, and by Q) and Q 
{B? + S 1>m + 2J lm bfBf ) )Bf ) ]nr n {R n 



p-i 



< 1 o H *l,m + $L,m X] & 1 — (i) ) 111 7 " ( 

i=l ^i 



<9ci 

2^(^(60-^) 



and 



(«/(£m - 5i, m )Bi 0) - 2J lm 0£ J ) lnr n (M m )| £ 
= -2J lm I J mm {n e b[ E) - «<JV(6i)) - n^y] 



-2J- lm a ( ? lnr n ( 



Substituted in (|6.6|), yields the identity 



gCo) 



(Bf + 5 lim + 2J lm 6^Bi u; )BW lnP n , F 4 



Wr(0)^(<) 
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+ \ (K e (e m - 5 1 , m )B? - 2J lm d£>) lnP n - Q, F ( 



si 



(o) 



This ends the proof of Proposition 6.1 



(6.9) 



(6.10) 
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For I = 1, . . . ,p, using the shorthand notation, 



Hp 



Fp 



g(0) 



-\Hf\F t 



(6.11) 



and ' := the equations (6.1) become (taking into account Yle=i ^ = ^) 



{X t ,Fi} = H t , l^E^p, with ^Xp = 0. 



£=1 



Proposition 6.2 Given for i = l,...,p functions Hp and Fp, such that the 
Wronskian of the derivatives F{, . . . , F' is non-zero, the system of ODE's 

{Xp,F e } = Hp, l^e^p 



subjected to the condition Yli=i — 0, has a unique solution (Xi 
Xp is given by 

i 

i 



, . . . , X p ), where 



( F[ F' 2 F' z 



Fp 

Xp = ^- det 
D 



F' \ 
p \ 



F'l F> 2 > F£ ■■■ -G 2 ■■■ F>> 



\ F[ p) F 2 {p) F 3 (p) 



-G r 



In this formula, D is the Wronskian of the functions F[, . . . , F' p , 

( F{ F 2 F3 

D := det 



F{' F% FX 



3 
1// 

2 x 3 



F' \ 

F" 

p 



(P) tt(p) j?(p) 



7^0 



1 A 2 - 1 3 

and i/ie Gi ? s are defined inductively as 

P tt 7-7(1) 



J 



Hp 



G i+ i — g\ + , c - 0, — ^ 



£=1 



Moreover 



det 



if* 


F 2 

F" 


^3 • 

F' 

r 3 

F" 

r 3 


■ F p 

. F' 

p 

. F" 

p 




r 2 


F (P) 
r 3 


F (p) 



0. 



(6.12) 



(6.13) 



(6.14) 
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Proof: If Xg is a solution of the equation {X e , Fi} = Hi, subjected to the 
condition Y^=i Xe = 0, then its derivatives are given by 



— Ggi + Xn - 



,(0 



(6.15) 



where for a fixed £, the Ge,i are defined inductively as 



Ge t o :— 0, G^\ :— — . . . , Gi >i+ i :— G' £ i H 1 —± 

re r n 



(0 



Indeed, starting with (6.15) and using X' e = y{Hi + X^Ff), one computes induc- 
tively 



X, 



(i+l) 



(j p z -+- y^ f —= — h — — — —5 — 



establishing (6.15). Summing up (6.15) for £ from 1 to p, one finds 



^=1 



Then solving this linear system for the X/s, one finds the ratio (6.12) above 



Then using that solution and expressing Y^t=i Xe = establishes (6.14) and thus 



the proof of Proposition 6.2 



B?\ with ' = #r = YZi Ju9 Ei , and 0£> with £Li 0£ 



This enables us to make the following statement, remembering the operators 



0. 



Theorem 6.3 ITie probability P n = P£(c,E) as in {l-Sj) , with the linear con- 
straint = ®> with Yl^=i K i = Ij satisfies a non-linear PDE in the 
boundary points of the subsets Ex, . . . , _E m and in i/ie target points \ . . . , 0^, 
given by the determinant of a (p + 1) x (p + 1) matrix 



det 





F 2 


F 3 .. 


■ F p 




r 2 


^3 •■ 


■ K 


if 


F" 

r 2 


F" 


■ K 




F <P) 


F (P) 
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G \ 



(6.16) 



G P J 



where the Fp, and Ce are given by in Proposition 6.1 and the Gg inductively 
by 



H 



(i) 



H 



(2)' 



Gn := 0. 



Proof of Theorem 1.1: It follows immediately from Theorem 6.3 by noticing that 



in the notation of (|1.11|), the Bf^ are expressed as 



(o) 



Proof of Corollary 1.2: The simplification comes from the fact that for one-time 



(i.e., m = 1) the operators do and dp differ by very little, namely: 



<9n 



d E , d e = dP + K f d E , e = e 



This means that the expression in brackets in the definition of Gj+i in (1.13) can 
be re-expressed as follows, 



rr(l) 



H, 



(2) 



° E - ^ ~f7 



Fp 



d, 



fr{l) rr(2) 

Fp 



(i) 



(i) 



KpHf, + 2KfbgF£, which one checks 



16 



to be the 



upon setting H i 

expression announced in (1.15) and one repeats the proof of Proposition 
with ' = 3e (instead of ' = do — —8e) and Xi \— > —Xp, ending the proo 



6.2 



of 



Corollary 1.2 



7 Examples 

7.1 One target point at the origin 

In this case, m = 2, p = 1 and the diagonal matrix A = 0. The matrix J reads 



J 



1 -c 
I - \ -c -1 



and one checks 



k = -1, bi = 0, «i = 1, d. 



(o) 



<9, 



(i) 



-n 9 9 



16 



Upon using the commutation relation [e E , d E ] = —d E and e — £q = 
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So, for £ = 1, one has 



TIC \ 

F 1 = -d d 1 \ogF n + — = -^Jd Ei +cd E2 ){cd Ei +d E2 )\ogF n + l 



nc 



(i) 



= 

and thus 



d 1 

-£ 2 <9 logP n = (e E2 - c ^)^3^ + c9 eJ l °S^n 

d 1 

-£ ^i logP n = -(e Ei - c -g^YZ^ ( c ^i + logP " 



(7.2) 



Go = 0, G?! = d 



leading to the PDE, with <9o and d\ as in (7.1) and H[ and Fj as in (7.2): (see 
and 0) 



7.2 Target points with some symmetry 

Consider non-intersecting Brownian motions leaving from and forced to p target 
points at time t = 1, with the only condition that the left-most and right-most 
target points are symmetric with respect to the origin, with p — 2 intermediate 
target points thrown in totally arbitrarily; this example will be used in section |8j 
It is convenient to rename the target points fix < ... < j3 p , as follows: 



a < —cx < . . . < — Cp-2 < 



(7.3) 



n p -2 



71- 



with the corresponding number of particles forced to those points at time t = 1. 
Using the change of variables ( 1.4[ ) from /3j's to 



b = (bx, b p ) = (a, -c x , —c 2 , -c p _ 2 , -a), 
one is led to the diagonal matrix of the forrrj^} 



(7.4) 



n+ 



ni 



rip-2 



A : = diag(a, . . . , a, -c a , . . . , -c 1( . . . , -c p _ 2 , . . . , -c p _ 2 , -a, . . . , -a) (7.5) 



17 Note the c, have nothing to do with the couplings appearing in (1.3 1 
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with the obvious constraint Yli K & 

1 



|a + \{—a) 



0, as in (1.5), and thus 
and Ki = for 2 ^ i ^ p — 1 . 



Moreover, setting c = (ci, . . . , c p _2), formulae (1.11 ) become 



5, 



(i) 



d_ 
da 



-d c ) ,d . 



j(p) 



2U 7' 6 



-, 2 < £ (7.6) 



and e = €e — a-§^ — £ c ; also set ' = &e- Besides the renaming n\ = n + , n % 



n. 



and rik 1— > n^^i for 2 ^ k ^ p — 1, already mentioned, one also has, referring to 
formulae (1.15), the following renaming: 

, F_, F k ^ F k _ u for 2 ^ A; ^ p - 1, 



+ 5 



ffM.-JZW Hjp^HW, h®~h®, hV>~h. 

1 — ► iZ^, #f 3 1 — ► i?"^, for 2 < £ < p - 1. 



(2) 
— ) 



Then, one checks from Corollary 1.2 formulae ( 1.15 )j_ that for 1 ^ I ^ p — 2 

8 



and for P := F*(E), with e = e E - a£ - e c (as in (|1.3|) for m = 1) 

1, 5 



2 -(T^-5 c + 9 E )5 B lnP+n ± , 



fi=^d B lnP+n A 



(i) 



d_ 
da 



-2d E e+(s + 3) -0 c + c> B MlnP + C ± , Hf> = 2 — InP+Q, 

^i 2) = ^(l-e±2a5 B )(T^ r -a c + o» E )lnP, iff > = (l- e -2c^) 5 - In P. (7.7) 



In accordance with formulae (6.11), adapted to the case m — 1, one defines for 
later use: 



{H { l\F ± } dB -{H?,F ± } x 

2 da "c) 



r(2) 



{H/ ,Fi} dls - {H^',Fe} 



(2) 



(7.8) 



and one checks that, with this notation (7.8) and upon decoding formula (1.15) 
for the Gfc's, 



Gk+i — OeGu + 



H+F_ 



(k) 



Fl 



+ 



H _ F (k) ^2 (fc) 



F 2 



F 2 



18 In the formulae below (7.7 1, the constants C± and have the value: 

/ 



C*4 



-n± I ±a ± 



P-2 

2E 



and Ci = 2ne 



p-2 



E 



V 



c^+a Q — a ^— ' C£ — c r 
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where is a shorthand for (dE) k F. With these expressions in mind, P := W$(E) 
satisfies the (near-Wronskian) PDE (1.12), i.e., 



det 



V 





F_ 




■ F p _2 


Go 


\ 




n 


F'_ 


F[ ■ 


F' 


Gi 






F» 


F'l 


F>> . 


F" 


G 2 




= 0, 




F (p) 




■ 1 p-2 


G p 


) 





(7.9) 



Special case: For Brow nian motions forced to a and —a, without the inter- 
mediate points, the formula (7.9) turns into the following determinant, with F± 



and H± as in (7.7), but with all c-partials removed: 



det 




= (H+F- 



H'_F + ){F + ,F_}=0. 



8 Pearcey process with inliers 



In this section, we consider non-intersecting Brownian motions leaving from 
and forced to p target points at time t — 1, with the only condition that the 
left-most and right-most target points are symmetric with respect to the origin, 
with p — 2 intermediate target points thrown in totally arbitrarily, exactly as in 
section 7.1. The purpose of this section is to identify the critical process obtained 
by letting n := n + = n_ — > oo and by rescaling a and the q accordingly, while 
keeping m, . . . , n p _2 fixed. This is the content of Theorem 1.3 



Proof of Theorem 1.3: The proof consists of letting n = n + = n_ — >• oo in the 
kernel (1.7) and in the PDE ( 1.12[ ). In the proof, which requires several steps, we 
shall restrict ourselves to m = 1 (one-time), except for Step 2, which deals with 
the kernel. 

Step 1: The PDE. The probability P := F^(E) satisfies the (near-Wronskian) 



PDE (7.9); see section 7.2. 



Step 2: The scaling limit of the Brownian kernel. Non-intersecting Brownian 
motions leaving from 0, such that n r particles are forced to (3 r at time t — 1, are 
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given by the kernel (1.7), which is, in this instance, conveniently rewritten as 



H tkl e (x, V\ 5, -Ci, • • • , -Cp_2, -a)dj/ 
dy 



, . dV dU - — - 

2n^(l-t k )(l-t e )J c Jr L U-V 



t f^- + PT L -n + \n{V-a)-n_\n{V+a) p-2 



e 
0. 



1-ti ~ r l-t» 



-n_(- ln(t/— a)— n_ ln(C/+a) 



r=l 



£/ + c r 
^ + c r 



dy 



y/ T{te—tk) 



for tk ^ 
for t k <ti. 



11) 



One then uses the same steepest descent method as for the case without inliers; 
the so-called steepest descent F-function is the one (depending on U or V) ap- 
pearing in the exponential, with three consecutive derivatives being = at the 
origin; the change of integration variables U = U'{n/2) 1 ^ and V = V '(n/2) 1 / 4 
then leads, in the limit for n = n + = n_ — > oo about the saddle point, to the 
kernel (1.18) (see for instance [T§\ and in the asymmetric case [3]). So, the limit 
is 



lim H nl V',a,ci,..., c p _ 2 , -a)dy 



h 
x - 

V z 



- I _|_ T k 
2 4V2n 
X 

4(n/^)V4 



Kl (X,Y; Ul ,...,u p _ 2 )dY, 



4(n/2)!/4 

a = \fn 



h = u e (f ) 



1/4 



where T .(X,Y; u±, . . . , Up-2) is the Pearcey kernel with inliers (1.18). 



1.2) 



Step 3: The scaling limit of the PDE. As mentioned, for the proof we limit 
ourselves to the one-time case, i.e., m = 1. We now proceed in two steps: 



(i) The change of variables (1.4) (especially footnote 6) from the non-intersecting 



Brownian motion probability to the matrix model (1.3); this change of variables 



appears in the first column of the table ( |8.3 ) below. In other terms, it is the time- 
dependent change from the variables ( ) to the variables ( ), yielding in 



particular the diagonal matrix A as in (7.5) 



(ii) Subsequently apply the scaling given by (8.2 ) with z := n 1//4 and a very small 



renaming s 



2 1 / 4 Wj, £ := X/2 5 / 4 for computational convenience. 



35 



This appears in the second column of table (8.3) below. 



x 
a 



X 

a 



t(i-t) 
it 



t 

x 
a 



4V2n 



X 



4(n/2)V4 



i.3) 



n 



Concatenating these two scalings leads to the following; in the string of equal- 
ities below, the change corresponding to (i) is indicated by =, whereas the second 
change (ii) is indicated by =: 



P(E,s,v) 



In Pi a '~ c 2>--->-c P -i, 

111 K n 



-a) 



(all Xi(t) E E) 



t= \(l + SZ 2 ) 

a = l/z 2 

£j = Vi/(s/2z) 

E = Ez/y/2 



In 



entries of diagonal matrix A\ 



E< 



t(i-ty 



2t 
l-t 



a, c, — a) 



entries of diagonal matrix A \ 
. " 

V2> 



2zE ll + SZ 2 fy/2 



^l-s 2 z A ' V l-sz 2 \z 2 ' z ' z 2 ) 



\ 



) 



: lnP~(£7; cvv-a 

entries of diagonal 
matrix A 



--: Q(E';a,c) 



t = l(l + sz 2 ) 
a = 1/z 2 
di = Vi/(\/2z) 
E = Ez/V2 



•4) 



Note that in the rest of this section, E and E' refer to complement of compact 



intervals; i.e., we shall be dealing with gap probabilities. The identity (8.4) 
suggests the ^-dependent map: 



T" 1 : (£, S ,^)^(£',a, Cj ), 
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given by 

E' = 

with inverse map 
given by 



2zE 



V2 /IT 



y/l — - s ~ : 



2^4 - 



1 — sz 2 '' 5 2 V 1 — sz 2 ' 



Vj 1 1 + sz 2 



i.5) 



T z : (-B',a,c) !-»• (E,s,Vj), 1 < j < p - 2, 



-£/ = ; ; r, S 



aV - 2 



a 2 z 4 + 2' " z 2 {a 2 z A + 2Y Vj az 
Then summarizing the above, one has 

Q(E';a,c) := logP^E'; a, c, -a) = logP^Tj 1 ^; s, =: P(E; s,v) 

and thus 



Vic- 



16) 



Q(E',a,c) = P 



( \[2azE' 



2^4 



a z 



V2~Ci 



U 2 z 4 + 2' z 2 (a 2 z i + 2Y az 

sati sfies the PDE ( 1.12 ) in the variables E', a, c, in terms of the operators specified 
in (7.6), with F±, Fg, H± \ given by (7.7). In order to express the PDE 
in terms of the function P(E; s,v), one must express all partials of Q(E'; a, c) in 
terms of partials of P(E; s,v) in E, s, v; e.g., 



d E ,Q(E';a,c) 



V2azE' 



d F P 



VI - s 2 z* 
2z 



d E P(E;s,v) 



a 2 z 4 + 2 

and thus the operators cL and d E , as acting on Q and P respectively, and similarly 
for the others, are related by the following; we also indicate what the relationship 
becomes for z — > 0: 



19 



<9 



,2^4 



V2 



d_ 
da 



s z a 
2z— d ° 



'l-sz 2 
z\l -d 



1 + sz 



2 "i 



z-sz 3 + ^s 2 z 5 + 0(z 7 ) ) (),, 



(1 - sz 




2\2 



1 + sz 2 d 
1 — sz 2 ds 







ds ~ Z2{£l ' + S ¥s ] - 8Z * 



1 SZ 2 . 2 , 

T+s~z^ {£v + SZ £e ' 



1 9 i - j 

S— + £ B - £„ + 0(2; 



2 9s 



J.7) 



19c 



Since 9 C = 2 9 Ci and 9„ = J2i Z ®vii the third relation is valid for <9 C and d v as well. 
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For notational simplicity, derivatives will often be abbreviated in the obvious 



way: 







{djFi^F^, (d E yp^p^, -P^P,... 



(8.8) 



while keeping in mind from (8.7) that cL acting on functions of (E',a,c), as 



F±, Hi and G?, translates, to leading order, into d E /(2z) acting on functions of 
(E,s,v); also notice the big gaps in the first few terms of the series for d ,. In 



view of the PDE (1.12), one needs the series expansion in z of the F's, the H's 



and the G's and their derivatives d E '- This is the content of: 

Lemma 8.1 Introducing the expressionY , withe = e E —e v , andv = (t>i, 

-Y := 4(e - 2s— - 2)d 2 E P + 16d v d E P + BP + {d E P, d 2 p} & , 



, Vp-2) 

(8.9) 



one checks, (remember from (7.8) the definition of H± and Hi) 



d E ,F e 



2z 



\d 2 E PT^d E P 



Be 
2z 



4y/2z 

>2 



0{z 



\d Vl d E P + n t - S -^-d Vl d E P ) • ' 'J'-' 



p-2 



H+ H— v — -\ Hp 

— — H h > — 

p ' p ' p„ 



64z 5 



p-2 



Hjd E iFi 



Y-3(d 2 E P)(d E P)) : Oil) 

(d%p) 



3 
32^ 



<9^ l+ ' 
2z 



P+0{z 



and also, for k = 0, 1, ... , one /ias 



G k+1 + 3 -^(F_-F + y k ^P" 



y( fc ) 



+ o(z~ fc - 1 ). 



16 



16(2^ 



,fc+2 



(8.10) 



(8.11) 



Proof: The formulae (8.10) are straightforward computations; one of them in- 



volves the expression Y introduced in (8.9). The big gaps in the series (8.7) of cL 



is responsible for the mere action of {d E /2z) 1 , in computing higher derivatives. 
Moreover, in the third formula, one notices that the sums 

hand side of = actually do not play any role in the leading terms, because Hi and 
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both are 0(1). Formula (8.11) is shown by induction; namely for k — 0, one 
checks, using the formulae (8.10), 

Gt + - F+)'P" 



16 

t,* + ^ F -- F +y p " 



16 



1 



6Az 2 

Assume inductively 
3^2 



(Y - 3P"P") 



3\/2 

16z* V4\/2 



pit pi 



O(l) 



Y 
64^ 



+ 0(1). 



Gi + ^-{F_-F + )^P" 



Y( i_1 ) 



16 



16(2«] 



— + 0(0 for 1 < z < fc, 



(8.12) 



and prove it for i — k + 1. Then, using the general definition (1.15) of Gk+i in 
terms of Gk, formula (8.12), the derivatives d^> of F± as in (8.10) and the last 
formula of (8.10), one checks 

G k+1 + 3 -^(F_-F + y k ^P" 



H_F1 




p - 2 H f Fl k) 



3V2 
~16~ 



(F_ - F + )( fe+1 )p" 



3 ^(K-F + )«P" + 0(^ i 



+ 



H + F 



(k) 



+ 



16 

H _ F (k) 



y( fc ) 



16(2*) 



fe+2 



0(z 



Fl 



^^^ + ^(F_-F + )( fc+1 )p" 

£=1 £ 



16 



establishing Lemma 8.1 



By Corollary |L2] Q(P';a,c) = lnP^(P';a,c) satisfies the PD E fl7.9 [), which 

and J8~5]. 



induces a PDE for P(E; s, v) = lnP„ (T z 1 (E; s, Vj)), remembering (8.4 



As pointed out, the PDE for Q(E';a } c) misses to be a Wronskian by the last 
column. It is appropriate to do some column operations; e.g., subtracting the 
first from the second and then adding the second, multiplied with P", to the 
last one; also it is convenient to multiply the columns with 2's and y ^'s. This 
gives us the determinant below, which vanishes according to Corollary L2) The 
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seco nd equality 
8. II above, 



uses in a straightforward way the series expansion of Lemma 



(F_ - F + )P" \ 

(F--F + yp" 

(F--F+)"P" 



= det 


f 2F + 
2F' + 

2F'l 


V2(F- 
V2(F_ 
V2(F. 


-F + ) 

-F+y 
-f+y 


2F l . 
2F{ 

2F'{ . 


■ 2F p -2 

• 2f;;_ 2 


G H 
GH 
G 2 H 


3^ 

16 
3^ 

16 
3V2 

16 




V 2F { ? ) 


V2(F_ 




2F[ P) . 




G P H 


3^ 
16 



(F_-F + )^P" / 



det 



p" 

(2*) 2 



(2i)3 + ^(^2 ) 



p' 

(2z) 

p" 

(2z) 2 
P'" 
(2*) 3 



O(z) 
-0(1) 



§£ + 2n 1 + 0(z 2 ) 

1 OP" 
(2z) 9«i 

i ap;" 

(2z) 2 Ot>i 



0(z) 
,0(1) 



V (2z)P+ 2 + ^IzP+l ; (2z)P+! + ^IzP-lJ ( 2z )P 8WJ +"^( z p-2j 



dP' 



1 



-27lp_ 2 + 0{Z 2 

+ 0(z) 



bp" 



(2z)dv p -2 



1 Qp(p+i) i V 

(2z)p 9t) p _ 2 ^zP-s J 



TeW + ^l) 



16(2z) 3 



y^- 1 ' +o(^ 

16(22)?+! T ^^P- 1 - 



\ 



c 

,6+p(p+l)/2 W P 



ds ' 



«9|P, . . . 







dv. 



p-2 



-d 2 E P,Y 



Ol 



1 



z A+p(p+l)/2 ■ 



The last equality = stems from the fact that the matrix consists of columns with 
increasing powers in 1/z, except for the element (1,1), whose leading term is 
2/z 4 . Therefore the leading contribution of the determinant of the matrix will be 
given by 

2 

— x the determinant of the (1,1)- minor, 

which indeed leads to equality =. Also the term —8sd E P n could be removed 
by adding 8s x (the first column); but we prefer not to do this, in view of the 
Taking the limit, when n 



conjecture 1.5 



oo, leads to the PDE for P 



limlogP ra . In the end, one must undo the slight renaming (|8.3|) of the variables 

v, = 2 l / 4 w 1 - 



T 



'.I- 



^i/2 5 / 4 an d go back to the (r, Uj, ^-va riables 



yielding §Y = 8 3/2 X, with X as defined in (1.20). This yields PDE (1.19), which 



ends the proof of Theorem 1.3 
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Very sketchy Proof of Corollary 1.4' A detailed proof appears in Adler-Orantin- 
van Moerbeke [3]. In the absence of inliers (p = 2), the Wronskian (1.19) is the 
determinant of a 2 x 2 matrix: 



= W 2 [d*d T \nF v , X] g = {d^d T la¥ p , X}. 



(8.13) 



Performing the same scaling limit on an asymmetric situation, with 2nq particles 
forced to — y/n and 2n(l — q) particles forced to y/n for < q < 1, with q ^ 1/2, 
leads to a PDE for the leading term having the form 



{$ \nF v , X} =0. 



(8.14) 



Thus lnP p satisfies two different PDE's, (8.13) and (8.14), given by two Wron- 
skians of X with d 2 d T lnP^ and d 3 lnP^. Then a functional-theoretical argument 
explained in [3] implies X = 0. ■ 



For inliers, we further conjecture -in analogy with the result in Corollary 1.4- 



the validity of equations (1.22) and (1.23), as stated in Conjecture 1.5 



9 Appendix: evaluation of the integral over the 
full range 



In this section we prove formula (|6.7|), i.e., we show that 



n 



)c=i fc=i 



1 1'/o, 



i=i 



k=l 



with g„(c) a function, depending on ci, . . . , c m _i and n only, computed below. In 
view of the representation of the above integral as the determinant of a moment 



matrix, as in (3.11), it suffices to prove that 



det 



(9.1) 



where, for 



M« : = 



p, the ri£ x N matrix JtfW is defined by 



\\dw k 



k=l 



5^ i < ri£ 
^ j < AT 
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Introducing for o,&6K the zero moment 20 
m(a,b) : = 



e -KEir=i u 'i- 2 Er=i c kWk w k+1 )+ awi +bw m j^dwk 

k=l 

(27t) m/2 - 



det J e -^ Jlia2+2Jlmab+Jmmb ^ 
we can express all the entries of as 



M. 



3* 5* 

da^ db* 



m(0,&f } ). 



(9.2) 



(9.3) 



Let us first prove (9.1) in t he cas e in which all ri£ are equal to 1 (so that p = N). 
Then, it follows from (9.2) and (9.3) that, for £ = the vector M™> is, 



modulo a constant which depends on ci, . . . , c m _i and N, but not on £, of the 
form 

Jf»>~ e -W (l,^^),...,^^)), 



where ay (6^) is a polynomial in b^p of degree j, with leading term ^— J7i m fef^ 
and whose coefficients are independent of £, but depend on c. It follows that, if 



all ni are equal to 1, therj 2 ^ 



det 



(9.4) 



proving (9.1) in that case. Let us show how the other extreme case, where there 
is only one ri£ (so that p = 1 and n\ = N), is derived from it. Let / : R — > R w 
be a smooth function and let /3 G R. Then 



det 



V/ (7vJl) (/?)y 



lim 



det 



/7(ft)\ 

/(A) 



(9.5) 



20 



Using 



e -4(Oi».™>+(^> dtUl ... du , = 



^/2 



for Q := - J- 1 and £ := (o,0, . . . , 0, b). 

21 Note g N (c) = (- J lm ) N{N ' 1),2 {- det J) N / 2 (2n) Nm / 2 , as easily follows from the argument, 

Tp T-rne — 1 



and finally in the full case g' N (c) — Jlfci Y\k=x ^-9n(c). 
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as follows by writing each f(fa) as a Taylor series around /(/?). Applied to 
/(/?) :=e-*P k f f, (l,a 1 (p),...,as- 1 (P)), 

> /5jv 



and 
then 



bf\..., we conclude using (9.3 ) and (9.4) that, when p — 1, 



■,(*) 



AT-l 



det (M (1) ) 



lim | [ fc! g N (c)e "a"* ^=x^ = g n (c)e 



(1)" 



ft,.. .,frv^f>r fc=i 



proving (9.1) in this case. The proof of formula (9.1) in the intermediate case, 
when there are several ri£, which are not equal to 1, follows in a similar way from 
(9.4), taking the limit fa — > b[ , for % = 1,...,N, with ni of the fa going to 



b^, namely fa, . . . ,f3 ni -> and /3 m+ i, . . . ,/3 ni+ „ 2 and so on, where 

now one divides by a product of p Vandermonde determinants, each going with 
a collapsing group. 



■(2) 
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